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ABSTRACT 

Based on sample data representing five years of 
monthly circulation totals from 50 academic libraries in Illinois, 
Iowa, Michigan, Minnesota, Missouri, and Ohio, a study was conducted 
to determine the most efficient smoothing forecasting methods for 
academic libraries. Smoothing forecasting methods were chosen because 
they have been characterized as easy to use and fairly accurate. It 
was found that smoothing forecasting nnethods worked very poorly on 
monthly library data due to the seasonality present in monthly 
library circulation totals. The only method recommended for use with 
monthly data was Winters' Linear and Seasonal Exponential Smoothing 
method, which has a specific seasonal component. Much greater success 
was achieved by using smoothing forecasting methods with 
yearly-lagged data, for example, using the circulation totals of past 
Januarys to predict the total of a future January. The One-Month 
Single Moving Average was found to be the most eflEicient smoothing 
method for forecasting future monthly circulation totals on 
yearly-lagged data with litle or no trend, while Brown's 
One-Parameter Linear Exponential method (with alpha set at 0.5) was 
recommended for use in trending yearly-lagged data. These methods 
ranked first and second respectively in minimizing both the mean 
percentage forecasting error and standard deviation of forecasting 
errors. A 27-item bibliography and plots showing the circulation data 
from the 50 libraries are included. (ESR) 
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Thxs :i>tudy idtutiixed forecasting metnoils as missxny from 
tne luahciyeaieiit sKilis of acdaemic liDrarians, Smootning 
forecdstiny utethoQh; were proposed as potential remedies. 
Furocat.tin'J researcm^rs nave characterized these tecnniques as 
easy tu use ana as lairly accurate for tne time and effort 
invet^ted m tneir calculation. They seem to De good candidate 
toiecastmq tecnnigues tor acaaemic liJ^rarians to try ^ A 
torecastinj competition was used to determine the most eificient 
smootninci xorecastmy methods for academic libraries. Tne 
saiuf).ie uatd were live years of monthly circulation totals from 
litty dcademic iiDiaries. Plots of each library's data snowed 
that tnt vast majority of libraries have heavily seasonal data. 

Tne results snowea that smootniny forecasting metnods work 
very poorly cn monthly library data ana we recommend that 
liDraiiahs not use tnem on monthly library data. The reason lor 
the poor performance of smoothing methods is the seasonality 
present in monthly library circulation totals. It a liorarian 
does wisii, however, to employ smoothing forecasting methods then 
we urge nim to use Winters' Linear and Seasonal Exponential 
Smoothing which nas a specific seasonal component . 

nucn l; rearer success was ach ieved by using smoothing 
forecasting metnous witn yearly-lagged data. An example of 

•■^ 3 



iacj'jiny aa^a one year is usxny tne cxrcuiation totals or past 
Januarys to predict tne total ot a future Januairy. It appears 
that lor lUciny libraries taKijiy a yearly-laygea approach 
completely de-seasouaiizes monthly circulation data. For some 
iioraries only a trend component is ielt in the data. 

tit: recDuimenaea that lor those libraries with little or no 
trend in their yeariy-iayged data/ the One-ttonth Single hoving 
Hverayc De used loi lorecasting tuture tnontnly circulation 
totdis. This method was the luost efficient sBioothmy method in 
minimiziny outn tne :uecin percentage error and standard deviation 
of foiecdstiny errois. For trending yearly-lay yea data, we 
rtf conaatfn ded that acaaemic librarians use Brown's One-Parameter 
Lmeai Exponential method (with alpha set at O.b) for 
JLorecastmy lature monthly circulation totals. This method can 
adjust its tortcasts when a trend is present. It ranged second 
m inmiruizmy ootn tne mean percentage forecasting error and 
standara deviation ol lorecasting errors. 
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fluny people have coutii i)u ted to tne ^Jiccess o± tnis 
enaeavGi:. Tjit^ iar9ei>t yroup ol unndmect ;i uniieraidea 
suppoiters dre tue iibrdriaiis who. took tj.me trom their busy 
scheuules Lo un^et our ciemands ior data about their libraries^; 
Tnost^ pc;op-Le nave oui many tuaiiKs. 

The Couacil on Library He:DOurces provided our tinancial 
suppor-: ana tnus lu a very real way made tnis study possible, 
both Dale i>entz^ Libiarian ot the University ol lowa^ and Cari 
Ot^jv^n, i/irector ot the bcnooi ot Library and Int ormation 
ScieJiCt:.* proviaea us with support and encourayement • 

Thi^ repc^ct could not have Deen written without the daily 
efforts of Judy butciier who spent long hours betore a computer 
terminal • . ■ 

i'inc.iiy our yratituae yoes to our families who generously 
forgave our enthusiasm xor library statistics through these many 
months. 
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INTHOjJUCTIuN 

There ijds ueen d rt^naissaiice oi interest in quantitative 
lurecaSLinq m science, management, and economics mtn only 
liaatea tranbter oi tnese new iorecasting techniques to the 
mandyement oi acacletaic libraries, forecasting is sucn a 
weii-(lt:veiopeu tool o± JDusiness that advice has been olfered to 
mandgers doout wnicn torecasting tecnnigue to choose [ Cnambers^ 
rtullicx, 6 limith, 1^71 J, Unfortunately, acaaeaic librarians 
don*t en^oy equivalent sage advice aoout wnicn torecastmg 
method to uiie. Tnis study proposed to fill this gap m the 
K'nowledye uase ol acaaeiaic librarians' by applying forecasting 
tecnniyues oi the smoothing variety [ MaKr idaKis 6 Wheelwright, 
1^78, cnapter J J to a large sample of acaueuiic libraries* 
Uiontnly circulations in a forecasting accuracy competition. The 
yodl Ol this research study was to tind the smoothing 
lorecaiitmg metnoas mat are tne most efticient forecasters ol 
acaaeruic lj.brary circulation data. / 

Li^^idry literat ure reveals little awareness ol t ort^castxng . 
insteaa, there die iaany statements made about the absence of 
quantitative tools in library decision maKing. Many authors 
lament tne current slate of library statistics. They have been 
chardcteriz^ed as primitive ny l^ogers and Vueuer [1^71, p. 275], 
and as DUsywurK by litrner [19b7, p.47J, burns [197^] has 
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cnaiTacte riz.ea tne cox-Lection oi iijordry statistics dis a routine 
operation with anciedi purpose that produces data that are 
utilized at a very unt>opnit>tica rea level. h6wever. Burns 
anticipated a i)etter future: "There will be a more 
sopais tica tea use ol statistics* to measure, lorecast, simulate, 
ana uioaei ail phases oi library operations especially tnose ot 
circulcttion" L 'l^V'/, p-/j. 

bviCKGhUUND UN FOiikCASl'I NG 

riambuiy [ 'iy7b, p-Jbj stated that "in library planniny and 
decision -maKiny predictions are mvariaiJly required". riis bold 
statement hc^s not Uiutivated iwucn theoreticax work or practical 
application oi turecastiny methodolocjies to library statistics, 
ii^ven as the t=conociic environment ol libraries has worsened, 
library jiianayurs nave not turned to tecnniyues such as 
forecastmy that would serve to tine tune a library 
act. .ties. "In particular, as resources dwindle libraries waut 
to predict tne Dehaviors oi their users and p'^tential users so 
tnat they may jjoth plan ana promote their activities" [Uervin, 
1^77 J. 

Tne use ot lorecastiny m library administration practice 
IS m sn arp contrast to tne acceptance ol forecasting m other 
aisciplinet>. iforecastmy , or trend analysis, is considered as 



an iiitejrax part oi bcxontitxc management and rationax decision 
liiaKxnvj. liaKCiaaKis una wneelwriyht descriDe torecastrng as a 
tool tiiat aeiiiiitci manayement to suioid an organization rrom tne^ 
vagaritri^ oi chance ev^jnts ana oecome more oietnodical in aeaiing 
wita Its eiiViLunuieiit J. li^Va, p.!:>j. LiKe oureauciacies 
everywnere^ acaaeiiiic iioraries neea tools that will ^nnance 
planniny and rational decision maKiny. One tool to nelp 
accoinplisli these manayerial tasKs is forecastiny. 

ForeCdStiny sauula De of interest to liJDrarians ana 
mioiiuation scientists xor at least two reasons. The iirst 
reason stems irorn tne aryument aLout tne desira JDili ty or 
manayeiiai rationality. LiKe ail manayers^ the library manager 
must allocate ais scarce resources prudently^ and make ''his . 
decisions based on his predictions of the effects of. alloca4:ing 
varyiny auiounts oi resources to the different functions in the 
library" [house^ l^'/^j. Prediction methods can potentially 
Decome one of the daily tools of tne library manayer. 

The second reason is more theoretical. Library -output 
statistics such as circulation data are intrinsically 
mterestiny variables that merit their own investiyation . 
i'revious wopk [brooics^ lyaij has demonstrated that 
library-output statistics nave surprising characteristics that 
are unanticipated i>y the rolKlore about them. i-'orecas tiny 
studies are only one methodology tor studymy library-output 
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stdtiij ticti • by demonstratiny an dDility to lorecast 

lii^rai >-output statistics^ we prove thatc^ we unaerstaad some o± 

tne oynainics tnat are arivmy them. 

LxDraiy literature is not distinguished by sophisticated 
applications ox torecastiny techniques . Many authors [Carpenter 
o Vasu^ ib^/d; noadiey 6 CiarK, l^Vz; Simpson^ i97bj writing 
aDout (juantitative or statistical methods m librarianship 
ijnore lorecastmg cumpietely. The topic is not treated in 
LancasteJL's The Measurement and Lvaluation ol Library Services 
and m anotiier volume^ Invt^ tiuative Methods m Library and 
Intgoictti oj' science: An Introduction , Wartyn and Lancaster- 
cover only the Delphi technique which xs^ method tor divininy a 
consensus 01 opinion. Conspicuously absent are any inier^utiai 
Statistical ioiecastiny ttchniyues. Rutherford Rogers and David 
Weber described the manageria 1 use. 01 liorary statistics as 
primitive and then proceeded to prove it by discussing 
torecastiuy only m terms of the descriptive method ot plotting 
trends on charts L*>^'/1# p.k7y]. They neglected to discuss any 
mtereiitial tecnniy ue that could establish the statistical 
signiiicance o± a graphic trend line. Stueart aad EastlicJc 
[lydlj^ who treated lorecastmy m three paragraphs^ also 
recommended only graphical methods. 

There are^ h<ftwever^ two other library lorecasting studies 
ot special note, liiriara DraKe [1976] considered linear 
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regre^iiiioii ajb d preaictxve tecni>icjue. She concluata that 
straight trena liuts are not the most efiicieut predxctors in 
dil library situations. Tne reason is that 'liurary data^ 
especially circulation data, show monthly or seasonal 
t luctudtions* Cyciicity may De one o± tne reasons that 
torecastiny techniques nave nad a retarded application to 
linrary statistics. Cyciicity of monthly library totals 
certainly played a .larye part in this study. The reader is 
invited to peruse tlie many monthly circulation plots given m 
Appendix J. rtost ot these data are strongly cyclical in that 
patterns that appear in one year are often repeated m other 
years. 

The most sophisticated forecasting 'study in library 
literature to date is oy Kang [1^79 J. He forecasted the 
requests for mterlibr ary loan sei vice^ received by the Illinois 
Research ana herei^ence Centers trom 1971 through 19?b using ^ 
several methods, including methods that can model cyclical data^ 
ana round regression to oe the best predictive technique. He 
used a weighted regression formula that gave less predictive 
value to older observatiojis, and greater weight to the most 
recent ones, liis study is flawed oy the tact that only one set 
or data wai:, used; hence, the generalizaoility of i^ang study is 
severely limited. 
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hi:;LAT£D F0K]:;CAS11NG STUDIES 

There are a jiuioDer ot rei ated stuaies that have attempted 
to predict cxrcaxation wxth causal techniques. Hodowanec [1960]. 
used multiple reyression analysis with twenty independent 
variables to predict circulation patterns ot graduate students^ 
undergraduate stuaeuts^ arid iaculty. McGrath [197b~ly77j 
isolated twelve ludepenaent variables in a multiple regression 
analysis to predict circulation of monographs by academic 
suDjects. Zweizig [197JJ used tne related approacn of multiple 
ditscriuiinan t analysis to isolate factors that determine public 
library use.. 

Another approach to predicting future library use is 
modeling standard statistical distributions on a sample of 
library circulations, for example^ LazoricK [1970 J found the 
demaua tor books in a collection to tollow a nt^gative Dinomiai 
distribution. NoziK [ 1974 j used a Karkov process to model book 
demand^ anu burreli [I9b0j offered a model to show likely 
patterns of future use of individual boox. titles. Morse and' 
Chen [ 197b J showed now Dias can De controlled in predicting the 
total yearly circulation of eacn class of books in a liDrary^ 
and Siote [1970] studied the past use of individual booKs as a 
predictor of tneir future use. Other investigators have used 
random samples to predict the total number of patrons entering a 
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libidry [Jones, ib^Vjj. 

It IS Clear irom tnis examination of related library 
literature tnat uo stuay nas ever used smootniny torecasting 



BeyiiiiJiuy in ueceiaDer 1^8/ arid continuing in January 1983, 
titty acaatiuic iibiaries from the American Library Directory 
were cnosen randomly from tne midwest states of Illinois, Iowa, 
hicnigan, Minnesota, Missouri, and Ohio. A copy ot the letter 
sent to eacn library is yiven m Appendix 1. 

Tne data reyuested were total monthly circulation counts 
tor u live-^year period, i.e., oO consecutive uiontnly circulation 
totals. The aim was to collect a set of time-series data from 
eacn liorary. Usaijle data were received from two libraries in 
Iowa, SIX libraritrs in Illinois, eight libraries in ttichiyan, 
one iit>rary m Minnesota, nine libraries in Missouri, and 
twenty-toui iiDraritrS m Ohio. 

Tnese data were loadeu in a computer ana a plot was drawn 
for each liDrary. Tne nionuily circulation plots appear in 
appendix The reaaer is invited to consult the plotted data. 

It IS our uDservation from tneso plots tnat m most instances 
academic libraries huve strongly seasonal patterns in their.. 



inetnoQs on acaaeiuic liorary circulations. 
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montnli^ cxrcuiatxon totals. That xs, it is evident in most of 
the piotb tnat certain characteristic patterns repeat tnenrselves 
trom yedL to year. 

SnOOThINC; fUKiiCASTING WliTHOLiOLOG Y 

Tne stuay applies the smoothing torecasting methods o± 
MaKriddKis and Wheeiwrigut ^ly/BJ to a sample or academic 
liDrary circulation statistics. These authors present formulas 
tor a numi)er oi smootniny methodologies. 

A siiiootninj lorecastmg method uses the intormation 
supplied oy previous data to create a torecast lor the tuture. 
It IS assumed in these methods that a sigjial exists in the past 
data that may ue oDscured oy a certain smaii percentage ot 
random errors, a smoothing torecastmg metnod weights certain 
past oust rvations and averages past random errors in order to 
reveal tne underlying signal in the data. Tnere are many 
metnods and many wayi> to utilize the information of past 
observations. It is often the case that eacn method permits 
many variations in either the selection of weights or in the 
time lags used tor smoothing. We have capitalized on tuis 
flexibility by systematically varying weights and time lags. 

bee a use oi tne uigniy seasonal nature of montnly library 
statistics, we employed the following strategy: all the methods 
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were run on tne moiituly data, and tnen trie metnods viere run on 
data laggtiQ one year, ihese y early-lagy ed data were developed 
irom tne 6u data points avaixaijle to us from each library and 
consistea. ot aatd points in positions 1, 13, ^b, 37, and 
Tnese aatd points were treated as a separate, ana smaller, time" 
series irom each iiDiary. -Not every library donated data that 
J:)egan wi tu January, thus txie lagged time series don't represent 
series thdt jump iroiu January to January to January, etc. Using 
data lagged m this rashion de-seasonalizes a monthly 
time-series Lecduse it transforms the data from a series ot 
months that run from January to December throughout one year to 
a series of the xdt:ntical months across several years. Thus it 
doesn't mdtter, lor example, that July may be a quieter month 
than August. With yearly-lagged data, Julys are compared to 
each other, as are Augusts compared to each other, etc. 

Following is a listing of the methods used in tnis study 
and the techniyues usea to initialize each method. 

St ylistic note: In the tollowing formulas braces have 
been usea to denotti subscripts. Thus f {t j should be read as f 
subscript t. 
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^ • Single Ploviug Averaqe ::^ 

Eicjiit varxationi:. or sin9ie movinq averages wer^L; used. An 
exposition oi tne sinyie movmy average tecnniyue is given by 
Wneelwriijht and ftaKriaaKis [lv7b^ P-^i>]- 

1. Une-wontii woviuy Avera^ye 
Ftt+lj^-Aitj 
Initialization: 

Note: Last montn»s total is used as next montn^s lorecast. 

. Iwo-nonta Moving Average 
1- It^lj = [X {tj -^X It-lj J/2 
liil tialization : ^ 

F Minx {ij t^j yz ^ 

Note: The average or the two preceding months is useu to 
forecast . 

3. Tnree-ilontn Moving Average 

F It 4 1} X it} +X It- Ij -i-X {t-2j J/3 
Ini titilization : 
Ft^i^LXnj^XUj^XlJj J/3 

Note: ihe average ot tne tnree preceding months is used to 
forecast. 
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4. rour-flonta Woviuy Average 

Fit^lJ^L^ltJ^Xit:-lj+Xtt~2j^Xlt-3j J/4 

Note: ihe dVeraye or tne roar precediny months is used to 
torecdst. 

5. Oiie/ lz.-'ioiitn Moviiivj Average 

F It + l^j ^X (tj 
liii txalizd tiou: 
Ftljj-Xllj 

Note: i^eriuittiiiy a month's total this year to oe the 
torecdst tor the sduie month next year removes the eltect of 
sedsondlity m the aatd. 

b. Two/ 12-Kionxh Woviny Average 
f It^^i^i-L^ Itj +x it-i/ij i/z 
lax tidlizdtxon: 

Note: The tirst and thirteenth observations are averaged 
to proviue a iort;Cast tor tne twenty-titth oDser vation • 

>• Tnree/li^i-honth Moving Average 
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F {t ^l^J =[ X Itj -t-X {t-rzi -i-X lt-24j J/3 
Initialization : 
FlJ7j-^[Xlzbj+X{lJJ+X llj J/3 

Note: uoservatxons xroin three precefe^xny years are 
dVuraged to proviae a iorecast tor the tourtn year. 

b. f^our/lz^-nontn noviny Averaye 

F (tt^lzj =[X itj -^X It-lzj -t^X lt"2Uj -^X It-Jbj j/4 
iai tializa tion : 

F 14 X iJ7j +X {zbj +X ll.Jj +X ilj j/i* 

Note: ODservations from lour preceeding years are averaged 
to provide a torecast tor the tifth year^ 

^ • i^J>ue a r novinq Average s 

Four variations oi the linear moving average method were 
used. Aii exposition or the linear moving average technigue is 
given ny wneelwright and MaK.ridaKis [ 197b , p.bbj 

1. Two-rtontii Linedi fiovin9 Aveidge-Lay One Month 
It + lj -a itj +D ttj m 
where m = 1 and: 

D Itj =Z lb' {tj -b" Iti) 

a It] -zS • (tj -S" Itj 
S"itJ=Lii' Itj+b' lt-1j J/2 
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b • It} ^[X ItJ +X It-lj J/2 
Initialization: 
S'l/j=LXUi-»-Xl1i J/2 
S- IJj =[X IJj +X l^i J/2 
S»lJj =[ (Jj +S' I2j J/2 
a i3j =^S • 13} -S" {3J 
b iij =^ (S' I3j -S" IJJ) 
F14) =alJJ 

Nott;: In tills aietnod two observations are averayed 
smootht-'U. This uv^raye itself is then averayed with a 
previousiji caXcuxattiu averaye to create a QouDle-siuootnea 
eftect. 

2- f n Ltjfc-i-lonth Linear ttovmy Averaye-Lay One nonth 
f {t + lj =u itj itj 
Where: 

i)ltj-iS' ttj -S"ltj) 

ci {tj - 2b' Itj -b" itj 

S"itj^[b' ttj It-lj +S' {t-2j J/3 
S • { tj =[ X Itj +X It-lj +X lt-2j J/3 
Initialization: 
b'l3j-[^lJi+Jil^)-»-XtlJ i/3 
S't4j=[Xl4j+XlJj*Xt2j y3 
b 'IbJ =IX[:>] +X I4j+X I3j J/3 

20 
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Nore: inih; metnod employs doubie sinootaxiig. Three sets 
three consecutive mouths are averaged, and the resultiucj three 
ctverctyts are aVeraoed to create a douDle-smoothed etiect. 

J. Four-i'ionth Lmeai woviny Average-Lay One nontn 

F lt-»-li Itj +0 ttj 

altj-^S« {tJ-S"ltj 

S"iti ^[b • it] +b« lt-1J +S» tt.-2j +S« lt-3] j/4 

S» (tj X It) -►X tt-1 J^X It-^i -►X {t"3j j/U 

Iiii tidiiz^ation : 

s» taj =[X i4j -^x {J J +x +x {1J J/^ 

o • ibj X ibj +X 14} ^jL lij -frX \2] J/4 

S • {OJ =L ^ tt>J +X lb] +X t^J -t-X {3} j/4 

5 • l7 J -[ X l7j -rX IbJ -i-X lb) -i^X {4} J/4 

S'H'i =[ S» 17} +S« lb} +S • ibj+S« 14} J/4 

al7}=2sM7j-S"l7} 

i>{7}=V^l^'l7}-b"l'}) 
Fl8}=a{7}-i^b{7}/ . 

Note: This metnoa smooths the tour preceeding monthly 
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observations to ciucite lour averages. These tour averages are 
dverayed tneiaseives. 

4. Two-Won tii Linear Moving Averaye-Lay Twelve rtontns 
F tt + l\cj =a [ti -^D itj 
bltJ-z(S»ltj-S"ltJ) 
a {tj lt>-S"itJ 
S"itj^Li5' ltj+S« tt-lzj J/-^ 
b'ttJ-L^ltj^Xlf-lzi J/2 
Initialization: 

S»t^bj=XXi2bJ +X {13] j/2 

a {kbj ^2::i » t2!:>j -"S^izSj 
D t^bj ~2 (b» t2bj i2b} ) 
F {J 7 J =a iz:>j +D {ZbJ 

Note: by usmy observations that are one year apart, 
seasoriaxity is eliminated. This metnod is identical to the 
Two-t^ionth Lintiar nov^ng Average above except that the uata used 
m tins metnoa are ragged one year. 
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C . Sang i i; Ex ponential bmoothing 

Eignteen variations oi the single exponential smoothing 
tecnnique were triea. Tne first nine varia tiont> are oased on 
shiitmcj the aipha wer-yht trom U.I to O.y by 0.1. The second 
nine variations incorporated the shitting alpha weiynt and 
employ data that are layyed one year. 

1. bmyle Exponential Smootniny (tne lirst nine variations) . 

F {ti- I J =F Itj +aipna (X [tj (tj ) 

wnere alpna steps i^y 0.1 Irom 0.1 to O.y. 
Initialization: 
Fl2j=Xl1i 

^M^ j ll] +alpha (Xt^j-Xtlj) 

i><ote: Tne aitierence l^etween tne last torecast and 
ooservation are yiven mrterent weights m detetmininy the next 
torecast. Tne yap oetween torecast and observation contributes 
the least wnen alpha is 0.1 and aiakes a larye contribution when 
alpna is O.y. 

2. single exponential Sifloothiny (the second nine variations) 

F {t->-1zj -Fltj -i-alpha (Xtt}-Ftt}) 

Where alpha steps £>y 0.1 Irora 0.1 to 0.9 
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lliitldliZatiOli: 

K l^bj ftj +axpnd (X tl^} ~X llj ) 

Note: ia*i data are layged one year to miaiinize tiie ettects 
Ox i>easondixty . 

U . br own's u ue-Paiameter Liiieai Lx ponential Smootumq 

IV "~ 

higiiteeSv Valid tioai^ ol brown »s. One-Parametiir Linear / 
Exponentxax binoothmy were tried. The tirst nine variations are 
baSeJ on sniitiny tne alpna weiynt irom 0.1 to U.y oy 0.1. me 
Second nine variations lucorpor a tea the shitting weight and 
employed datu that are lagged one year- 

1. brown's uue-Paraiueter Linear Exponential Smoothing 
(Lne lirst uiiie Variations) 

F Itf Ij =d it] {tj m 

Where a =- 1 
i> itj -=^Lcilpiia/ (l-alpna) J (S • {t} -S" {t} ) 

wnere aipna steps oy 0.1 trom 0.1 to u.y 
aiti=^S» ltj-b"ttj 

Si» (tj -alpna (tj (1-alpna) S" [t-'lj 
b»{tj=aipha X itj t (l-alpha) b • tt-1} 
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S'{1j-Xi1} , 

by2j=:.xpnci X l^j + (1-dlpha) X ilj 
S'H^i =^dlprm 5 • l^j + (l-dlpiid) X 11 j 

-aipna/ci-aipud) iZj -6'' [z^ ) J 

Nott: brown's urie-Parameter Linedr Exponential Siaootning 
niethoa IS a composite ot the Linear Moving Average method ana 
tne bingle Lxponentidl Smootning method. As is tirue with tne 
exponential i;>UiOOthiny methods, heavier weights will maKe this 
metnoa respond aiore yuicKl^ to rapidly lluctuating seasonal 
data. 

2. rrown^'ij one-Pa la meter Linear Lxponential Smoothing 
(the secona nine variatons) 

F It+l^j =altj+/{t}m 

Where m 1 

W 

b lti=Lalpna/(l-dlpiid) j(S' lti-S"{tj) 

wheie aiphd t>teps by U.l Iroiu 0.1 to U.y 
a it J It] -5" l i-i 

b ' {tj ^axpud S • IL] ♦ (l-aipha) S" [^-^2] 
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S»itj ^dlpnd X (tj + (1-aj.pna) S» {t-lzj 

Initidiizdtioii: 

S'l1j=X{li 

S"il]=Xili 

S ' 11 j) =aipnci X ilJ>+ (1-alphd) X {1} 
S"l13j=aipha S« tl33 + {l-alpha)X{l3 
a i13J =ZS • {I3j -6" tUj 

b{lJj-^Lalpha/(1-alpha) ](S« lld}-S"tldj) 
F 12bj =a {13] +D lldj 

Nott: iivery twexith observation is used to minimize the 
ertect oi seasonality. 



^- brown's ^uaaratic Kxpopeutiai Smoothing >. 

tiynteen variations ot brown's Quhdratic Exponential 
baioothiny were testea. The rirst nine variations are Dased on 
shiftmvj tne alpha weight ±rom 0.1 to 0 by 0.1. Tne second 
nine variations incorporated the shitting weight and employed 
data that are lagyea one year. * 

1. brown 'hi v^uadratic Kxpooential Smoothing 
(the tirst nine variations) 
F tt+1j =CL Itj {tj ffl + 1/2 c {t j (m squared) 
where m - 1 
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c {t j dlpiid syudrea/ ( l-alpn€i) sguarea J 
(S» Itj -kbi" (t) -^S* • • [t] ) 

wnere cilpiia stepped by 0.1 from 0.1 to 0.9 
b ItJ -[dipha/k (1-alpha) S4uared J 

L lt>-balpna) bi* (tj - (lU-balpha) S" (tj 

(4--jalpha) s • • r{tj ] 
a ttj =JS « Itj -jS" itj +s • • • ttj 
S • ttj =diph<i S" tt-i ♦ (1-aipha) S' " [t-IJ 

(tj -alpha S • ttj + ll-dipha) S" tt-1} 
S • ttj -dlpiia X ttj + ( l-alpna) S • (t-lj 
Ini tializa'tioii : 

bMij-xtlj 
s"tij=xtij 

s"»tiJ=xtij 

S'tzj=alpha X t^j ♦ (1-aipha) X tij 
S"t2j=alpna S • t^j ♦ (1 -alpha) Xtlj 
S • ' • t^) =aipiia S" t-ij •»■ (l-alpha) X 1-1} 
a t/j " l^j -iS" tzj +S" {2j 
b t^j =L axpha/^ (l-aipha) squared j 

[ (b-bdlpha) S* tzj - (10-8alpha) S" {2} + 

(4-3aipha) S • " {2} 
ct2j =Lalpna scjuared/ ( l-alpha) squared J 

(S* t2j-^S"t2j+S" '{2}) 
F t3 j =-a t2} +0 t2j + \/2 c t2j 
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Note: brown's yuadrcitic ibethocl^ which is an extension ot 
linear exponential smootning^ locuses on trends that are more 
complex than imeai trends. Triple exponential smoothiny is 
used* AS witn previous weiyntiny methods, heavier weights ailow 
this uiethod to tracK rapidly fluctuating data. 

^. brown's v>uadiatic Exponential Smoothing 
(the second nine variations) 

i {t -Hzj =a ttj ni ^ b [tj 1/^c {t) (in squared) 

where m = 1 
c {t} -[ alpna squared/ ( l-alpna) squared J 

(:d» {tj -2S"tt.)-»^S» • • ttj ) 
b it j =[alpna/^ (1 - alpha) squared ] 

[ (o-balpna) S • {tj - (lU-tialpha) S" {t J + 
- ti4-Jaipha) S • » » (tj J 
a it J ^ • it} -3S" tt^j -i-S • • • (t J 
S • • • iti =^alptia {t j (l-alpha) • • lt-12J 
S'»{tj=alpha S • {tj ♦ (1-alpha) S"{t-12} 
S»it}=aipha X [tj + ( l-alpha) S • tt-1i2j 
Initialization: 

S"«i'iJ=Xi1j 



EKLC 



2S 



23 



S»llJj=aXpha X {IJj ■^ (I'-alphd) X (1j 
S"l 1 Jj ^alpiid tlJj (l^cilpha) X {1j 
S • • • t "^-^J =^alpiict b" t '^•^i ^ I l-alpha) X t Ij 
a 11 Jj 113] -JiS" tlij +:5 • • » {13} 

^ -^i -[ aipna/-6 1 1-aipha) squared ] 

[ ib-balpha) S» fI3j— ( 10-balpha) ll-^j + 

It4-j5aipha) b « • » ll3j J 
ctl3j=^[alpha syuarea/ (1 -alpha) squared J 

{b« ll3j-^b«» t13J+S»» • tl^i ) 
F 12 bj ^a ll3J-+b t13J +l/2c 1 13} 

Nottr: r-vcry twelftn oJDservation was usea to minimize the 
effects of seasonality. 

• iioit^s Twu-i^arameter Linear Exponential Smoothinq 

Uiyhteen variations oi this method were tried. The first 
nine Varied aipna ana gamma Dy O.u irom 0.1 to 0.^. The two 
parameters tooK the values 0.1^ 0.5^ and U.9 consecutively 
resultm^j in nine comoinations . The second nine variations 
repeated this approach with data lagged one year. 

I. holt's Two-Parameter Linear Exponential Smoothing 
(the tirst lime variations) 
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F (t -fr^lj ItJ Ui 

wnere lu = 1 
b (t j ^cjciuiiaci (5 It j tt- Ij ) ( 1"9dmiiia) b {t-lj 
S [t j =aXpna X {tj (l-alpftci) (i>[t-1}+ b {t-lj 
Illl tiaiizdLxon: 

bl1J=L (X(2j-X[1})^(X{4j-X(JJ) ]/2 
S {2j =ajLphd X [2> t- (1-alphci) 

:X{1]+ ((Xl2j-X[li)^ (X {4j-X{J}))/2] 
o t2j =yaaiijia* (Si2j-Xtlj) + 

( l-gaiiima) [ ( (X I2j -X { 1] ) - (X {4} «X I3j ) ) /2 ] 
F iJj=Si2i+bl2j 

Note; jiolt'ia metnod iS ditterent trom preceeding metliocls 
because ol the necessity ot specifying two parameters: alpha 
ana yaaiaia. 



2. Holt's Two-Pardaieter Linear Exponential Smoothing 
(tne secona nine variations) 

F It-t- 12J It] +Dltj m 
Where m =^ 1 

b itj -yamiaa (S itj lt-12j ) + ( l-^amma) b {t-12) 
Slt}-alpna X [tj (l-alpna) (Slt-12} btt-12j 
initialization: 
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(AllJJ -xnj)-i"(X{J?J -X 125}) 3/2 
b {l^j =^aiphci X 11 J} + 

( l-alphd) L X llj 1 (X { 1 JJ -X { 1} ) + (X { J?J -X I2bj ) ) /2 J 
b 11 JJ ^9diaiuci (S (1 jj -X {Ij ) + 

(l-yaaiaia) L ((XllJJ-X (IJ) + (X 137} -X I2b}) )/2J 
F I2bj lljj +D [1 JJ 

iNote: Tne etiects of sed:ionaiity are controllea using 
data laygea one year. 

G . Aa aprive-tvesponse-fvat e Single Exponential Smoothip g 

iNine vaiiations of tnis metnod were tried. Beta was variea 
xroni 0 .1 to 0 . y Dy u . 1 . - ' — . .. ^ ^* • 



F Iti-lj =^alp]ia ttjXitj + ( 1-alpha {tj ) F It} 

waere alpaa (tj = I ii^lt} /M j^t} | 
fc;{t} = oeta e{t} ^ ( 1 -bet a) E lt-1} 
M ttj ^Deta I e{tj I + (l-Deta) W it-lj 
ett} - Xttj - F (tj 

initialization: 
fl2|=X{l} 
, e {2J -Xl2} -X{1} 
« l2j=i>eta|X(2}-Xll} I 




h l^j =-betd (X t^j-X llj ) 
d j.pha {2} -Detd 

i''lJj=L>eta Xi-^j -f (1-lDeta) X llj 

e tJj = X{j>-Ft3j 

tt IJJ =beta|XiJj I 

£ {J j =Dera (X {3 j {JJ ) 

alpha {Jj -Dera 

F j =DetaX l^j + F 13} 

iNote: This uietnod reyuires only the specitication o± the 
r>eta value. This method is adaptive in the sense that the alpha 
value will Change when therj^ is a change in the Dasic pattern ot 
tne aa ta • 

d • Winters ^ Linear and Season a l Exponent ial Saoothing 
Winters' method worKs on the development of three 
eguations. riach eyuation lipcuses on one aspect ot a pattern in 
time--ser les data: the stationary levels the trendy and the 
seasonality ot tne data. Tne values ot the parameters were 
syjitiuetricaiiy varied ±rom 0.1 to 0.7 by O.K That is, the 
parameters were assigned the same value so that^ tor example^ 
aipna yamiiia = Deta - 0.1. Later ^ other variations were tried 

in au eiiort to achieve lower errors. TaDle II shows the 

■ ■ • \ . 
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specilicaLxon o± tne paramtiter i» used. 

Jr' it + nij = (b (tj + X)(tjM)I {t-L+m} 

wfiere L, the lag, is one year 
It + lj ^ + ultj ) I it-12J 

S tt J ^ (alpha {X It j /± (t-Lj ) ) + ( ( 1-aipna) (s {t-1 j +d {t-l) ) ) 

0 ttj^ijamiaa (b Itj -Stt-lJ) + ( l-yamma) £) 

1 {tj = (Dcta (X Itj /s ItJ ) J + (1-beta) I{t-L} 

Initidlizatioii: 
S {1 Jj =X (IJj 

Dt13j=[ (Xl Uj -X ilj) + (X{14J-X{2)) + (X{lbj -X{JJ ) J/Jb 
1{1)-X{1J/X 

wnere X=suui A [1 to 'I3j/li 
F{l4j-(X11JJ t1Jj)l{2i 

^ • A Typical Trend J:;quation 

Tat:; tollowiiiy metnod was used as a general purpose 

equation loi trend. 
Flt-Hj = ^X{tJ-Xtt-1i . ' 
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l«ote: This is> act d smootninq eyuation described by 
fSdKridaKis diid Wne^riwriynt Dut is based on our oDservatioii o± ' 
tne Leducibility oi soiootxiing methods into several generic 
equdtiohi^. bee iippenuix Z for an exposition ol the reauciDilxty 
ot suioothmvj lorecastm^ methods. 

KESULTS 

Tne results oi tins study are presented in Table 1 through 
Table IV . 

A nalysis ot Tanies 1 > I I# and III 

The lirst three taoles are headed by columns lor 1. average 
percentage error, z. standard deviation (SD) , J. coetlicient of 
variation (CV), and M. a miniaization of the mean and standard 
deviation, nereatter Cdlled average torys number (AFN) • 

Ihe average percentage error is the percent error in 
forecasting atter tne maicated number of runs averaged for all 
bO libraries. In other words, all the libraries aie treated as 
It they were iiity Vciriations of one library and the fdrecasting 
results were averaged lor all fifty libraries. A measure of 
variability tor this mean is tne standard deviation (SD) . The 
coefficient of variation (CV) is the ratio of standard deviation 
to the mean. The CV provides a way of comparing the 
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characteristics ot the airterent aiethods. The avercige lorys 
nuniDer (ArN) is dti attempt to provide in one numoer an index o± 
the i^ize oi ooth the iuean error and standard deviation. A good 
lorecasting metnod wouid minimize both ot tiiese. Ihe 
Pytnagorean theorem provides a method ±oc doing this., The 
Pytnayorean theorem states tnat the square on tue hypotenuse ot 
a right triangle eyuals the sum ot the syuares on the other 
sidt^s. The application or tnis tneorem to create tne AJPN is 
based on our foliowic^ij observation. If the average percent 
error were grapncd on one axis and tne corresponding standard 
deviation were yraphea on the other axis^ each forecasting 
metnod could r>e represented uy a plotted point. The single best 
method would oe represented by that point closest to the origin 
ot the axes. Tne distance any point is from the origin or the 
axes IS the length ol the hypothenuse on the graph created by 
the aistances irom the origin oi the mean and standard 
deviation . The AFw permits tne iorecasting metnods to oe 
compared. Tne most eriicient iorecasting method would minimize 
JDotn Its average percentage error and its standard deviation and 
thus have the smallest AJb'N. 

Table I presents tne results ot running the smoothing 
forecasting methods on the monthly data provided by the fifty 
iioraries. The results are reported alter the torty -seventh run 
throug h the data so as to facilitate compar ison among tne 
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p^crtnods (x^e.^ rtintclii • inetnod aDsoroed U datd point '* iior 
initialization) . * 

bcverai ine tnods are distinguished in this y roup with low 
AFN'i;,. They are:' One-Month iloving Average^ all ox the Single 
exponential bmootniny Methods^ and the Adaptive-Kesponse-Rate 
Single hxponentiai Smoothing (JDeta - O.yy) .' As a group the 
average Ai-iN tor the rorecasting methods is 4b. -^y. It would 
appear tnat siiupier iorecastmg methods have the greater success 
witn monthly data. A real assessment^ however^ must JDe made 
against methods that can handle the seasonality ol library data 
such as Winters' method. " 

Table 11 presents the results ot runningf Winters* metnod on 
the mont]:ly aai:a supplied Dy the bU libraries. The average AfN 
for Winters • metnod is which is much smaller than the 

average AFN lor the smoothing rorecasting methods on monthly 
data. In tact^ tr.e Qiiierence m mean AFU is statistically 
signiticant (t = Z4.17^ p = U.UOUO). It is clear that winters^ 
metnod minimizes AFN with Diontniy data. Winters • metnod takes 
seasonality into account and therefore outperforms the other 
smoothing methods with monthly data. 

TaJjle 111 presents the results ot running the smoothing 
Iorecastmg methods on the yearly-lagged data lor the t>0 
libraries. The AFl^'s seem very much smaller. The average AFN 
tor Table III is 10. bJ. This is significantly smaller than the 
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average kFN lor Taijie I = 21.1b, p = O.UUUO) and also 
siynii icantiy smaller tnan tne averaye aFN tor TaJ^le II (t = 
J .iU^ p = O.UUib) . 

II is clear rnat binootning lorecasriny metnods perlorm mucn 
bettui on yeariy-iay ytiu Ucita. The reason is ttxdt the 
y edrly --layyccl data itiinoves tne seasonality from the time series. 
The ae-seasonaxizea data doesn't stay at the same levex eacn 
year out remains t^teady enouyn tor tne smoothing methods to do a 
mucu better tracKiny job tnan tne smoothiny methods do on 
stronyiy seasonal data. 

The Two-nonth Linear Woviny AVeraye has the smallest aFN 
and uidv De reyaraea as the nest smoothing torecasting method on. 
t he average . It iuust be emphasized that in Tables 1 to 111 that 
the results are averaged across the bO libraries tor each run. 
No one particular iinrary can expect to have or dia have 3.27 
AFN. l^^iis AFN IS an average lor the tifty libraries. 

It IS premature to conclude t^at the Two-Wonth Linear 
Movmy Average snouia i^e recommended tor all libraries. There 
IS a need to evaluate the many tluctuations ot the mean percent 
error ana standard deviation values that occurred during the 
many iorecasts made tor each method. Before we recommend any 
torecasting inetnoa^ we want to investigate the behavior ot each 
method tnrouynout all its torecasts. For an appreciation ot the 
Vctriability ot forecdsts other calculations are needed that 
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accummuiciLe tor ecast errors throuyh tne forecastiny cyc±e tor 
edcn metnoa. Tnis vrouiG give a i^etter picture ot tne ability of 
the metnoa^ to peiiotm with indiriauajL libraries. 

A nalysis oi laPle IV 

TaUie IV presents a comparison ol tne forecasting m^tnods 
with y early-iayyea aata. New comparison calculations haVe been 
used: 1. tne dveraye ttean Syuared hrror {i1S£) , 2. tne average 
Mean Absolute Percentaye Error (ilAPE) , 3. the average Mean 
Percentage Error {s^tt.) , *4 . tne average Standard Deviation (SD) , 
and b. tne Average Forys Number (AFN) . 

MdKriaaKis and wneelwriyht detine these averages in the 
loliowing way (p.bbb-8t^) : 
M ean Sgu ar eu Error ( MSE) 

The mean scjuarea error is a measure or accuracy computed by 
squaring tne mdiviauai error lor eacn item in a data set and 
then rmumg tne average or taean value ol the sum ot ol those 
squares. Tne mean squared error gives greater weight to large 
errors tnan to small errors Decause the errors are squared 
betore being summea. 

fj ean A us o iute Percentage Error (M APEj 

The mean absolute percentage error is tne mean or average 
ol tne sum ot all ot the percentage errors tor a given data set 
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taictfii without reydia to siyu . (Ihat is, then absolute values 
are summed aiiU tiie civeiaye computed.) It is one measure ot 
accura cy couiiiiojiiy usea m quantitative methods ot toiecasting. 

M ean t^cicenttigt: tttor (fIPE ) 

Tne luear* percentage error is tne average ot all ol tne 
percentage errors ioi a given data set. This average allows 
positive ahd negative percentage errors to cancel one another, 
because oi tnis, it is sometimes used as a measure o± bias m 
tne application ot a lorecasting method. 

TdOle IV IS unique in giving iniormation about the 
variauiiity m iorecas ting errors ior each lorecasting raetnod* 
Tne three averages delinod above were cumulated throughout the 
t c ..ocast liiy series lur each library lor each method and then au 
average was touna. The average ol these averages across lifty 
libraries is presented m Table IV. This table consequently 
presents average errors, standard aeviatibns, and AFN's that a 
liDrarian may expect to Una it ne were to use any ot these 
methoQs to forecast his own library's circulation data. 

An analysis oi Taoie IV leads to the tollowing ranKing ot 
smootning lorecasting laetnods with yearly-lagged data: 
1 ) Ohe~flon tn bingie no via g ii verage 

This ranKs tirst because it had the smallest kFH, the 
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smdixe::.r averayo hb^, tne iiecoud smallest WAPE, and une secoad 
smallest Wi-h. It we yuantify tnis pt-riorniance and give a score 
ot onu tux: aaviny tnti smalietit mean, a two tor haviay the second 
smallest mean, etc., tnen this method ''has a comparative ranX 
score orl-^l-rz-^-^^b. 

2) i:»rown*s one-Par aaieter Linear Exponential Smoothing (alpha - 
O.b) ' 

Tnis aietiioa ranKs second because it had tne second smallest 
average rJSK, the second smallest average AfN, and -the third 
smallest average riAPE and tourtn smallest C1P£. This metnod has 
a comparative ranK score ot 2 2 3 ^ ^ - ^^ . 

3) Single J^xponential Smootning (alpha - 0.99) 

This method ranxs third because it had the smallest average 
i^APE, tne third smallest average AFW and MSE^ cind the sixth 
smallest average fiPE. Thus this method nas a comparative rank 
score ot 1 + 3 + 3-^6 = 13. 

^) urown's quadratic Exponential (alpha = 0.3) 

This methoa ranks Jtourth ;>ecause it had the fourth largest 
rtSE, MAPE and AFi^ anu the lifth largest MPE. Thus its 
comparative ranK score is 4 + 4 4 5 - 17. 

b) Holt's Two-Pa. .uecer Linear Exponential Smoothing (a = 0.9, g 
= 0. b) ' 

Tni::i metjiua ranV.s iiftn because it scored tittn on average 
h^H, USE, ana dAPE . It lanKed third with rtPE. Thus its 

<\0 ' 
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comparatiVt ranK score isb+5 + 5 + 3= 18. 
b ) Two -Month Linear novxmj Average 

This mernoa ranks csxxtn as it had tne seventh largest tiSii, 
MAPE, AFN but the smallest nvE . Thus its comparative ranx score 
is 7 ^ 7 + 7 + 1 = z/, 

7) Ada pt ive-hcsponse-Kdte Single Exponential Smoothing (beta - 
0 .yy) 

This aiethuii v.nKs last because it had the sixth highest 
mean scores iOL . . the averages calculated except KPE rfhere it 
ranKea icist» Thus its comparative ranK score is b + o + b + ? = 
25. 
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TABLE I 

hLSULTS Of USING MONTHLY DATA 

Average % 
Error after 







4 7 Kuns 




cv 


A PN^ 

AX M ^ 


Single rtoving 


Averages 










One-nontn 




lb .04 






Oh 


Two-ttonth 




20 .Ub 


38 .b3 


1 92 




Tnree--ttorAtn 




22^77 


3b . lo 






Four-rt on th 




25.13 


34.98 


1.39 


43.07 


Linear Moving 


Averages 










Two-rtoiitn 




11 .bb 


bb. Id 


4 7 ■^ 


Sh 40 


TJiree-Month 




14 .07 


b2.4fa 


-J « « o 


SU H 1 


Four-Month 




lb .bl 


4b .bb 




S () A 0 


Single Exponentidi Smootning 








alpna -0.1 




22 .y7 


31.11 


1 hS 


•ift k7 


alpha = 0.^ 




23. b1 


31.b3 


■ . O / 


J? . O 1 


a 1 pn d = 0 • J 




22. b7 


32.4^ 






alpha = 0.4 




21.72 


33.05 




AU ss 


alpha = O.b 




20 .b2 


33.7b 


1 •b4 


3y .bfa 


alpna - o.p 




ly .bi 


34.4b 


1.77 


39 .59 


alpha 0.7 




18 .42 


3b. 02 


1 *yo 


39.57 


alpha - O.b 




17.37 


3b. 4b 


2.04 


39.4b 


alpha - 0.^ 




lb .36 


3 b. 71 


2.1b 


39.29 


alpha - O.t^y 




15. b2 


3b, 83 


2.31 


39.0b 


Brown one-- Parameter 


Linear Exponential 


">othing 




alpiia - 0.1 




Z4.27 


33. bb 


1.34 


41.41 


alpha = 0.^ 




22 .b4 


33.9 b 


1 .50 


40.81 


alpha - O.J 




20 .20 


3b .22 


1.7^ 


41 .47 


alpna - 0.4 




17.50 


39.0b 


2.23 


42.79 


alpha = 0 




15 .04 


41. bb 


2.77 


44 .28 


alpna = O.b 




12 .b7 


43. bl 


3.38 


4b.37 


alpna = 0.7 




10 .y3 


44. bO 


4.08 


4b. 92 


alpha - u.b 




9.22 


4b. 2b 




4b. 19 


alpha ~ O.y 




7.b5 


4b. 9b 


6.01 


4b. b9 


alpha - 0.^^ 




b .28 


47".U2 


7.48 


47 .44 



*An explanation of these abJtireviations used in Tables I 
through 111 is cuven m the iiesults section. 
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'MbLE 1 (Con«t.) 






urowii ■ 


s quadratic Exponential Smoothing 






alpna 




zJ.bl 


33. 5j 


1.42 


4 1.01 


alpna 


= u . / 


20 .bb 


3b. 7b 


1.7b 


42,16 


alpna 


= U . J 


lb . IJ 


42. oy 


2.bl 


45.07 


alpna 




^2 .27 


4 7.13 


3.84 


4b .70 


alpna 


= O.b 


y .yv 


4y.0b 


4.92 


bO .Ob 


dlpna 


= U .0 


0.75 


^'2.25 


7.74 


5 2 . b .'S 


alpha 


- u . / 


4 .bb 


54.0*1 


11.79 


b4 .20 


alpna 


= u.o 


2 .bb 


57.49 


21.44 


b7.bb 


alpna 




o.7b 


b3. 8y 


81.42 


b3 .b9 


alpha 




-1 .2b 


72. bo 


^-b7.bb 


72 .67 


Holt 's 


Two-Paranieter Linear 


Exponential Smoothing 




alpna 


= 0.1 yamma = 0.1 \^ 


bb .z3 


bb.by 


1.0 1 


94 . 13 


alpna 


= 0.1 gamoid - 0 .b 


J2.74 


4b. 18 


1.38 


5b. 80 


alpha 


= 0.1 y amma = U •'i . 


2b .21 


39.21 


1.39 


4 b .30 


alpna 


= o.b gam ma =- 0.1 


27 .bb 


37. J4 


1.35 


46.47 


alpha 


= O.b gamma r ^-b 


18.11 


41. bb 


2.30 


45 ay 


alpna 


= O.b gamma = 0.^ 


lb. 41 


51.01 


3.31 


53 


aipaa 


=• 0 .^j gamma - o. 1 


ly . J4 


^B.oy 


1.97 


*T ^ . # ^ 


alpha 


0 gamma =- O.b 


11 .82 


44.20 


3.74 


4 b . 75 


alpha 


- 0.^ gaiuma - 0..^ 


8 .27 


47.2b 


b.7l 


U 7 


alpna 


- O.^U gamma =• O.yy 


b .38 


47.13 


7.39 


U 7 Sh 


Aaaptive-Kersponse-i<a te Single Exponential 


Smoothing 




£>etd - 


0 . 1 


27. by 


3:^.bb 


1.1b 




Deta = 


0 


2b .35 


33.71 


1.33 


42.1b 


beta = 


0 .J 


22.17 


34.80 


1.57 


4 1 .2b 


Xitita = 


0 .4 


21 .lb 


3b. 37 


1.72 


42. U7 


Deta = 


O.b 


20.0b 


37. 8y 


1.8y 


42.87 


beta - 


O.b. 


18 .41 


38. b7 


2.10 


42. b3 


beta - 


0.7 


17 .Ob 


38.3d 


2.25 


42.00 


Deta = 


O.b 


lb .ly 


37.88 


2.34 


41. ly 


beta - 




15.81 


37.1b 


2.35 


4U. 3b 


beta = 


0 .yy 


lb.8i* 


3b. 9b 


2.2.7 


3y.2y 


Trena 


Lq Uation 


b .13 


47. 18 


7. by 


4 7.bB 
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TABLE II 



WlKTJiK'i) METHOD ON 


MONTHLY 


DATA 






Average % 










hrrors Atter 










47 Runs 


SD 


CV 


AFH 


Coefticients =• U.l* 


t).u2 


1U.30 


1.72 


11. yb 


Q.2 


0 .22 


11 .24 


1 . 8 1 , 


12.8b 


U.3 


b .82 


10.4:^ 


i.7y 


11 .94. 


U.4 


1 .ye 


14 .b2 


7.3b 


14 .bb 


U.b. 


3.01 


12.1b 


3.3b 


12. b7 


U .0 


3 . 7b 


11 .b4 


3.1b 


12.42 


0.7 


4 . 37 


12.28 


2.8 1 


13.04 


U.7 U.b U.?-** 


4.bb 


ll.bb 


-<;.4y 


1^.48 


u.y 


4.b2 


17.31 


3.7b 


17.92 



♦Each oi the three parcimeters tooK. this value 
(i.e., aipna = U.l, beta = U.l, yammd =^ U.l) 

**The tniee parameteis took these three values 
(I.e., aiuna = u.7, beta = U.b, gamraa = U.b) 
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DATA 
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Siiiyie rtovmg Averdges , 










une-nonth 




b.47 


7 


.40 






lwo--Mont h 


J 


V.bd 


b 


.oy 




9 • 9 b 


Three~«oii tli 




12.7b 


7 


.10 


0 .5b 


lU 59 


Two-ftonth Linear Moving Average 














2.33 


2 


.29 






biJigie Expoiientidi Smootiiiug 










dlpnd ' 0.1 




11. ZO 


y 


. j5 


0 .83 


14 .59 


alpne. =• u.z 


3 


10. bU 


b 


. b4 


0 . b ' 


13 80 


aipiia - 0.^ 


3 


10.01 


8 


.3b 


0 .8 ^ 


1 OU 


alpha - O.U 


J 




7 


.88 


0 .83 


12 .30 


alpha = O.b 


J 


ti.y 1 


7 


.46 


0.84 


11.62 


alpna - u.b 


3 




7 


. 10 


0 .84 


11-01 


alpha ^ 0.7 




7.t^8 


b 


• 83 


0.8b 


10 . bU 


alpha = u.c 


3 


7.6 1 


b 


• 64 


0 .87 


lU . 10 


alpha - 0.!:^ 


3 


7.32 


b 


.57 


0 .yo 


y . b4 


alphd = 0.!:^!:^ 


3 


7.13 


b 


.01 


0 .93 


y .72 


brown's One- 


Parameter 












alphd -0.1 


3 


10. oO 


b 


.84 


0.83 


13.80 


alpna = o 


3 


^. m 


7 


.84 


0.83 


12 .2i> 


alpna - U.i 


3 


d.27 


b 


.93 


0.84 


10.7'J 


alpha =• 0.4 


\^ 


7.23 


b 


.lb 


0.85 


y .bO 


alpha - 0 .b 


3 


b.33 


5 


.o2 


0.89 


b .4b 


alpna =- o.o 


3 


b.o2 


b 


.3b 


0.9b 


7.77 


alpha - 0.7 


3 


• b. lb 


b 


.42 


1.0b 


7.aB 


alpha = u .b 


3 


4.yb 


b 


.86 


1.18 


7 .bb 


alpna - 0.^ 


3 


b.Oy 


b 


.b9 


1.31 


b.Ul 
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TAbLE 


111 (Con • t .) 






Brown 


yuadratic 


Exponential 


Smoothing 






dipha - 


0.1 


j\ 


10.00 


a. 34 


0 .63 


u . u z 


dipna - 


0.2 








0 .b4 


1 . / J 


dipna - 


0. 


j' 


0 . 0 2^ 


b.bg 


w . o 


o . # o 


aipud ^ 


0.4 




b.27 


a. 93 




/ . 


alpha - 


0.5 




U.30 


4.7b 


1.11 


K ill 
O.HI 


dipild = 


0.0 




3.BJ 


b.22 


1.36 


^^ 11*7 


diplld = 


0.7 


J 




0.37 


l.bl 


/ .QU 


aipna = 


0.0 




4.82 


b.3? 


1.74 


? . O D 


dlpnd =- 


o.y 




b.bl 


11. Ub 


1.7fa 


1-3 . 1 # 


holt's Two-Parmetfer Linear Exponential Smoothing 




d=U, 1 


0.1 




m.bo 


13.40 


0.92 


■ -J . # H 


a=0. 1 g^ 


0 .5 




13.!^b 


12.77 


U .91 


18.92 
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0.^ 


J 


13.43 


12.1b 


0 .90 


18.11 


d=U.t> g = 


0. 1 




10. 3b* 


9.14 


O.titi 


13.84 


a^U.b g= 


O.b 


i 


d. 2b 


b.87 


0 .b3 


10.74 


a=^U,b g= 


0 .y 




0.3b 


b.0 4 


0 .79 


8.11 


d=0.^ g- 


0.1 




7.^0 


7.30 


0 .92 


10 . 7b 


a~U.b^ g= 


0 .b 


3 


b.ob 


b.bl 


U.99 


7,9b 


a = U.y g=^ 


o.y 




4.bU 


b.ob 


1.22 


7 .32 


a-0.!:^y g 


^o.gy 




b.Ub 


7.67 


1.U1 


9.41 


Trend ii^quation J 


b.60 


7.yy 


1 .43 


9.7( 
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TABLE IV 

^CunPAHlSUN UF XEARLi -LAGGED DATA rtETHODS 
iNumner /ivt-raye Avercye Average Average Average 
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41 


.99 


Thref-fiohth 


I 
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27 .Ob 


12 .7b 


59 


.00 .. 


bO 


.b3 


four -Hon tn 


1 


b n^.bi 


33 .3b 


22. b1 


87 


.2b 


90 
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Two rlonth Linear 


Movmy Average 
















3y74 .00 
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57 


.17 


b7 
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alpiici=0. 1 
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i4bd.67 


34 .69 


11. zo 


67 


.b9 


o8 


.6b 


aiphd=0.2 


S 
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31 .02 


10. oO 


61 


.13 


b2 


.22 
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J 
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bb 


.3b 


bo 
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dipild = 0. 1 


3 


2y2b.9b 


25 .bV 


9.44 


50 


.25 


b1 


.33 


aiphd=U . b 




24lb.b 1 


23.99 


8.91 


4b 


.84 


46 


.93 


aipna=^u. 6 




^1039. 4 b 


22 .4b 


b.42 


42 


.18 


43 


.2b 


dipUd=U. 7 


J 


1774.04 
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7.98 


39 


.28 


40 


.3b 
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J 


IbO 1 .by 


20 .3b 


7.61 


37 


.lb 


3b 


.22 


aipna-u.y 


J 


1b07.04 


19.73 


7.32 


3b 


.8b 


36 


.89 


dlpnd-U. 


3 
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19. bb 


7. 13 


3b 


.39 


J6 
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brown •s One-Parametei Linear Exponential Smoothing 



dipild=U. 1 


J 




3U 


.87 


10.60 


61.04 


62.13 


alpha=U. /. 


J 




2b 


.bl 


9.41 


bO.OO 


bl .08 


alpha=U. ^ 


6 




11 


.31 


b.27 


41. b4 


42.89 


alpna=0. 4 


:J 




20 


.80 


7.23 


37.00 


37.9b 


aipha=U. b 


J 


1448. yo 


20 


.4b 


6.33 


3b. b3 


30 .33 


alpna=U. b 


6 


lb3d.7b 


21 


.lb 


b.62 


36.72 


37.36 


aipna=^U. V 


6 


1734.26 


23 


.b4 


b.lb 


3y .48 


40.09 


alpha=0. d 


:J 


I^y3.b7 


25 


.b4 


4.^6 


42^98 


43. b2 


alpha=0. ^ 


3 


>:3b7.14 


27 


.32 


b.09 


47.21 


46.00 



*Tnese aDi)revid tions are explained in tht; Results section. 
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TAbLE IV (Con't.) 
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4 1 
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J 
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b.51 
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Exponential Saioothmy 
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38 .8b 
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87.80 


89 
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8b 


.4b 
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J 
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31 .34 


10.39 
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62 


.48 


a=U.b y=U.b 




^73b.b2 


28 .57 


8.25 


50.12 


51 


.05 


a=0.b y-o.y 




2z40 .b4 


2b .81 


b. Jb 


45.75 


4b 


.38 


a=0.^ y=0. 1 




2 J3b .4b 


2o .71 


7.90 


44.72 


4b 


.71 


a=u.^ y=u.b 




iyy4.47 


.48 


b.bb 


41.72 


42 


.40 


ci = U.*^ y^O.y 




z434 .bo 


28 .02 


4 .b4 


47.09 


48 


.23 


a=u.yy y=0. 
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/ 9 .tib 
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b3 
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20 


.b8 
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b9 
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20 
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b3 


.b4 


bb 
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U .3 1 
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29 .bl 


19 


.51 


bO 


.89 


b3 


.94 
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0.4 1 


3599.20 


28 .54 


18 


.bO 


57 


.04 


59 


.99 


beta 




U .5 1 


3'!7o.b5 


27 .^1 


17 


.bb 


53 


.48 


5b 


.38 


Deta 




U .0 1 


2b 19.02 


z .' ..,i9 


17 


.38 


50 


.17 


53 


.09 


beta 




0 .7 1 


^r+7.u8 


■ •.. .42 


lb 


.84 


47 


.58 


" 50 


.47 


Detd 




0 .tt 1 


2 37o „ l 


.09 


lb 


.39 


45 


.91 


4d 


.75 


beta 




U .9 1 


£. ^. '-/ i . Jh 


/.b . 19 


lb 


.34 




.07 


47 


.94 


beta 




U .99 1 


^2b2o9d 


2o .20 


lb 


.32 


4 4 


.b8 


47 


.57 


IT en a 
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6 


jutil.bl 


iO .4/ 




.bO 


5j 


.79 


54 


.79 



is 



CONCLUSION 

Tnit> stuay nas attemptea to assess smoothing torecasting 
methous as eilicient ^predictors ot academic library statistics. 
Alter fexaiuminy the results of this study^ we would liKe to maKe 
tne lollowing poiJits : 

1) we dou^t recoiumeua tne application ol smootniny torecasting 
metnods on monthly library circulation totals. The data^ as 
exniDiteu in the piors ot Appendix are simply tar too 
strdsonal lor siuootning methods to model them accurately . This 
conclusion is m accord with other assessments oi these methods^ 

1. e.^ naKridaxis and Wneelwrignt (1^7b^ P-t)b^) mh^o comment on 
inaDiiity ol smoothing iaetnods to handle seasonal data. As 
revealed in TaDle 11^ a me.tnud liKe Winters* wnich takes 
Seasonality into account does a better lorecasting job. but we 
recommend none or tnese methods ior monthly data. However^ it a 
\\uv rxi.\.xii were to employ smoothing forecasting methods then we 
would LiL^e mm to use Winters* Linear and Seasonal Exponential 
Smootiiiug. This method^ however, is relatively complicated and 
copjuuies axot ol data m initialization. 

2. ) we urge the use or smoothing torecasting methods on 
yearly ^lagyea data. As nad Deen snown. Table Ill's AFN scores 
are signiixcantly smaller tnan either Tables I ana 11 Ai?N*s. 
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But DC fort* deciain^ whicn ^mootuinq torecastiny metnod to use 
witn y eariy-xayyeu aata^ we investigated the variability oi 
errois that a iiDrariaa miynt expect to find with nis library's 
data. Taole IV reveaieu tr^at the uae^Month tloviny Averaye and 
Brown One-jearameter Linear Exponential Smoothing (alpna = 0.5) 
were tne two D^st methods to use on yeariy-iay yed library 
circulation totals. 

3) We uiye acddemic librarians to regard a plot ol monthly 
circulation data layyed one year and maKe the following 
decision: if their uata don't show trend across several years, 
then use tne Oae-nonth floving Average Uiethod for predicting 
tuture totals. li , on the other hanu, their data ao trend up or 
down across. i>everai years, then we urge the use of Brown's 
One-Pa ru meter Linear Lxponential Smoothing (alpha =• O.bj . 

The implications of these lecomiaendations are that with 
trendiess or nearly trendless yearly-rlagged data, the best 
predictor ol a uonth next year is the that month's total this 
year, witn treaaing yearly-lagged data, the best predictor is a 
aethod tnat smooths several previous yearly-laggea observations 
together as a trend aujustment. 

our reeling is thaL more research needs to be done on the 
pre-analysis stage or rorecastmg. By examining their aata 
before forecasting, academic librarians may ne able to select 




ttie most dppropLicite tore cds tiny met noil. We Oelxeve that, 
xmprovtia torecdsts wxii be the result. 
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liitr Letter ol Iiiyuiry 



In SepteniDer we Deyan a research project whose axm is to 
lind tne most etiective method ot predicting tuture circulation 
levels oi dCddemic liuraries JDased on past circulation data. We 
are writiiiy to you as part ol our collection pnase; we are 
looj<.ing xor academic lioraries that would have approximately 
five years' wortn or montnly circulation counts availaJ^le lor 
analysis iii our study. 

Our ijiteiit IS to collect data trom aDout 5U academic 
libraries in the ilia west. Each library would contribute five 
years' wortn of monthly total circulation counts (i.e., oO 
consecutive montniy total circulation counts) . Each time series 
tnus collected irom eacn iiDiary would tnen De analyzed with an 
interactive iorecasting software package called SlBYL/HUi4i^£R 
that is uvuiiaole ror research use on a Hewlett PacKard 20U0 
computer here at the University oi Iowa. From the output 
produced by SlbYi,/hUNN£K we will be aDle to determine which of 
24 extrnpolative time-series methods would be able to model each 
library's uata and muKe tne most effective forecasts. Our study 
of dcadeuiic library circulation statistics follows in the 
traaicion or other more general xorecastmg studies, e.g* "The 




Accuidcy ui LxLrcipoiatioa (Time StrKbs) riethoas; Results of a 
l-'orecctibtiiivj competition" uy Spyros WaKridaKis, et al. Journal 
ot Forec asting , v.l, payes lll-lb3, 1V^b2. This study has 
pionteiea ex trapolative techniques in yenerai; we Joelieve we 
will oe the lir^jt to apply them to academic library circulation 
counts . 

Can you i^ena us auout live years^ worth oi your toral 
montniy circulation countsi* Any consecutive sequence ol 60 
iiionths m tne recent past will do. We Know that liDraries 
collect uiiierincj statistics and call them ditterent names*. 
These uiittrtnces will not attect our study tor we aon"t intend 
to compare the counxt: irom library to liurary. Instead we want 
to study the perioraiance ot the z4 lorecasting methods on Uiany 
ditteront sets ol circulation data. 

An example ot tiie type ol data We are looKinq lor would be 
tne annual statisticdl summary many lit)raries compile yiving the 
total circulation ot the main and any branch libraries "Tor each 
month ot tne preceediny year. A pnotocopy ot sucn an annual 
statistical summary wouxd suit our purposes very well. 

Tile contribution ol your library's total montniy counts are 
important to our study and we would like to thanK you now lor 
every ellort you expend on our behalf. 
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APPENDIX 2 

Tiifc i<t:auciDility oi bmoothiny Forecdstiny Metnods 

Sutoothiny turecds ting methods torni a lamily ot methods 
that, when coeiticients are set to 1^ reduce down to two simple 
formulas: 

The Firs t Formula: 
F it^mj ~ A {tj 

dua When m = 1 
F tt + lj =^ X {rj 

This IS equivalent t-o sayiny thut a luture observation in a 
tiflie series will be liKe its immeaiate predecessor. 
T he beco nd Formula: 

F {t-^mj - X (t j + (X {tj - Xit-li)m 
and When m - 1 

F it 1 1 j - X tt J (X iti - X It- 134 

or 

F{t + 1j = ^X{tj - X{t-1} 

This IS ecjuivdi^-nt to saying that a luture oDservation in a 
time Series will be liKe it.s immeaiate predecessor plus the 
ditterence observed m tne time series between the last two 
observations • 

AS is illastratea nelow^ smoothiny torecastmg methods 
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degeiit/iate into either ot these two toriuulas wnen coetticients 
are set to 1 . 

S l ucj j.e hxpon ential Smoothing 

f (t-Mj - dipiia X {t j (1-dlphd) F{tj 

wnen dlpiia - 1 . 
F tt+ij =^ X (tj 

brown's une-jr^ar ameter Linear Lxb o nential Smootnxnc) 
S'ltj ^ alpha A Itj (l-alpna) S»lt-1j 
wfien aipnci - 1 
' b'ttj - Xitj 

b"it) = alpha 5 • {tj ( l-aipna), It-lj 

When aipna - 1 
S"tti - s« Itj 

Thereiore b" [tj = ^'Itj ^ X{tJ 
a ttj = jLY. itj X Itj 
a ttj = Xitj 

u Itj - \/\) U {tj - X itj ) 

oltj = U . 
F (t +mj =^ a Itj 

Wnen m - 1 and since a {t j ^ X {tj 
Flt^lj - A Itj 
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br own's (.udaiatic i:.xpoiiearial SiiiOothxng 

wheii dxphd f 1 

waicn red xxc^s^ to 
a (tj = A tt} 

Bit) - [ alpna/z I l-alpna) syuarea J 

[ (D-b axpha) S»ttj - (10-8 aipna) S"t ^ 
(M-J alpha) b» * » (tj J 
when drpna 'H 

i) ttj = o 

c ttj = Ldipiici squared/ ( 1-alfna) squdrea J 
■ V-^' It} -zi>"{tj + Itj) 

Wueii A±[)iia = 1 

i-' {t + mj = X itj 
When iu = 1 

i'tt+lj^Attj 
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ri Qlt*s T wo-pardinutei Linear i:;x poiifcatial Smoothiag 

S^tj = dlpna Altj (1-alpna) (S it-lj D(t:-lj) 
wjiea dipha - 1 

L){tj = gamma (:> {tj -Slt-lj) ^ (l-gamma) D{t.-1} 

whtr jj y amiua = '1 
b ttj - Sit) - b It-lj 

or 

D It] = X It] ^ X tt-lj 
F (t tn.J = b {tj ♦ i> ItJ m 

w n e li ih - I 
Fit-* IJ ^ itj - X {t-1j 

Aa ap ti ve -h tjs pun sfc- Kate Single Exponential Smoothing 
fit+1j = dipna ttj X{tj ^ (1-alphdttj) r [tj 
alpha (t^lj ^ liltj/Mltjl 
E (t j = Deta e it] ( 1-Deta) E (t-lj 
H It j ^ beta leltjl + (1 -t>eta) (t-lj 
G ttj = X ttj - r ttj 

wnen Deta - 1 
L{ti = elt] 
ItJ ^ leltj I 

tlieret ore 
alpna {t^lj =1 
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W xiiteri:>' Linear and Seasonal lixPQnentiai Siiioothxng 

S itj = dipna (X itj/I {t-Lj.) + (l-alpna) (S ix-lj + blt-lj) 

When aipiid - 1 
b{tj - XttJ/ilt-Lj 

b (tj = tjcimifla (S Itj - Sit-1}) + (1-yaniUia) btt-lj) 

wneii yamDici - 1 
b It) ^ {tj - S {t-lj 

I It] =^ beta (X ttj/S it] ) + (l-beta)'! {t-L} 

wnen beta - 1 
I icj = A [tj/i>lti 

initialize b{tj = X {tJ 
i {tj - 1 

b {tJ - X |t] - X It+lj 
f Iti^ Ij =/:X irj - A it-lj 



64 



by 



Circuxdtioa Plots ot i^xtzy LiJDraries 
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